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Abstract. Fix strictly increasing right continuous functions with left limits 
Wi : R -* R, i = l,...,d, and let W(x) = E^Li^C^i) for x e K d . We 
construct the VF-Sobolev spaces, which consist of functions / having weak 
generalized gradients Vyy/ = (&Wif> dw d f)- Several properties, that are 
analogous to classical results on Sobolev spaces, are obtained. W-gencralized 
elliptic and parabolic equations are also established, along with results on ex- 
istence and uniqueness of weak solutions of such equations. Homogenization 
results of suitable random operators are investigated. Finally, as an applica- 
tion of all the theory developed, we prove a hydrodynamic limit for gradient 
processes with conductances (induced by W) in random environments. 



1. Introduction 

The space of functions that admit differentiation in a weak sense has been widely- 
studied in the mathematical literature. The usage of such spaces provides a wide 
application to the theory of partial differential equations (PDE), and to many other 
areas of pure and applied mathematics. These spaces have become associated with 
the name of the late Russian mathematician S. L. Sobolev, although their origins 
predate his major contributions to their development in the late 1930s. In theory 
of PDEs, the idea of Sobolev space allows one to introduce the notion of weak 
solutions whose existence, uniqueness, regularities, and well-posedness are based 
on tools of functional analysis. 

In classical theory of PDEs, two important classes of equations are: elliptic and 
parabolic PDEs. They are second-order PDEs, with some constraints (coerciveness) 
in the higher-order terms. The elliptic equations typically model the flow of some 
chemical quantity within some region, whereas the parabolic equations model the 
time evolution of such quantities. Consider the following particular classes of elliptic 
and parabolic equations: 



(i) 



for t G (0, T] and x £ D, where D is some suitable domain, and g is a function. 
Sobolev spaces are the natural environment to treat equations like (fT]) - an elegant 
exposition of this fact can be found in [2]. 
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Consider the following generalization of the above equations: 




d t u(t, x) = Y,i=i d Xl d Wi u(t, x) 
u(0,x) = g(x), 



where dw t stands for the generalized derivative operator, where, for each i, Wi is 
a one-dimensional strictly increasing (not necessarily continuous) function. Note 
that if Wi(xi) = Xi, we obtain the equations in ([I]). This notion of generalized 
derivative has been studied by several authors in the literature, see for instance, 
[U HQl HU H2] ■ We also call attention to p] since it provides a detailed study 
of such notion. The equations in ((2]) have the same physical interpretation as the 
equations in JTJ). However, the latter covers more general situations. For instance, 
[Hj and |16j argue that these equations may be used to model a diffusion of particles 
within a region with membranes induced by the discontinuities of the functions Wi . 
Unfortunately, the standard Sobolev spaces are not suitable for being used as the 
space of weak solutions of equations in the form of (|2|) . 

One of our goals in this work is to define and obtain some properties of a space, 
which we call W-Sobolev space. This space lets us formalize a notion of weak 
generalized derivative in such a way that, if a function is VF-differentiable in the 
strong sense, it will also be differentiable in the weak sense, with their derivatives 
coinciding. Moreover, the VF-Sobolev space will coincide with the standard Sobolev 
space if Wj(xj) = Xi for all i. With this in mind, we will be able to define weak 
solutions of equations in We will prove that there exist weak solutions for such 
equations, and also, for some cases, the uniqueness of such weak solutions. Some 
analogous to classical results of Sobolev spaces are obtained, such as Poincare's 
inequality and Rellich-Kondrachov's compactness theorem. 

Besides the treatment of elliptic and parabolic equations in terms of these W- 
Sobolev spaces, we are also interested in studying Homogenization and Hydrody- 
namic Limits. The study of homogenization is motivated by several applications 
in mechanics, physics, chemistry and engineering. For example, when one studies 
the thermal or electric conductivity in heterogeneous materials, the macroscopic 
properties of crystals or the structure of polymers, are typically described in terms 
of linear or non-linear PDEs for medium with periodic or quasi-periodic structure, 
or, more generally, stochastic. 

We will consider stochastic homogenization. In the stochastic context, several 
works on homogenization of operators with random coefficients have been published 
(see, for instance, |13[ I14j and references therein). In homogenization theory, only 
the stationarity of such random field is used. The notion of stationary random field 
is formulated in such a manner that it covers many objects of non-probabilistic 
nature, e.g., operators with periodic or quasi-periodic coefficients. 

The focus of our approach is to study the asymptotic behavior of effective coeffi- 
cients for a family of random difference schemes, whose coefficients can be obtained 
by the discretization of random high-contrast lattice structures. In this sense, we 
want to extend the theory of homogenization of random operators developed in 
|14j . as well as to prove its main Theorem (Theorem 2.16) to the context in which 
we have weak generalized derivatives. 

Lastly, as an application of all the theory developed for IF-Sobolcv spaces, elliptic 
operators, parabolic equations and homogenization, we prove a hydrodynamic limit 
for gradient processes with conductances in random environments. Hydrodynamic 
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limit for gradient processes with conductances have been obtained in [6] for the 
one-dimensional setup and in [16j for the d-dimensional setup. However, with the 
tools developed in our present article, the proof of the hydrodynamic limit on a 
more general setup (in random environments) turns out to be simpler and much 
more natural. Furthermore, the proof of this hydrodynamic limit also provides an 
existence theorem for the generalized parabolic equations such as the one in l[2]) . 

The hydrodynamic limit allows one to obtain a description of the thermodynamic 
characteristics (e.g., temperature, density, pressure, etc.) of infinite systems assum- 
ing that the underlying dynamics is stochastic and follows the statistical mechanics 
approach introduced by Boltzmann. More precisely, it allows one to deduce the 
macroscopic behavior of the system from the microscopic interaction among par- 
ticles. We will consider a microscopic dynamics consisting of random walks on 
the lattice submitted to some local interaction, the so-called interacting particle 
systems introduced by Spitzer [IS], see also [S]. Therefore, this approach justifies 
rigorously a method often used by physicists to establish the partial differential 
equations that describe the evolution of the thermodynamic characteristics of a 
fluid, and thus, the existence of weak solutions of such PDEs can be viewed as one 
of the goals of the hydrodynamic limit. 

The random environment we considered is governed by the coefficients of the 
discrete formulation of the model (the process on the lattice). It is possible to ob- 
tain other formulations of random environments, for instance, in [3] they proved a 
hydrodynamic limit for a gradient process with conductances in a random environ- 
ment whose randomness consists of the random choice of the conductances. The 
hydrodynamic limit for a gradient process without conductances on the random 
environment we are considering was proved in [TJ. We would like to mention that 
in [4] a process evolving on a percolation cluster (a lattice with some bonds removed 
randomly) was considered and the resulting process turned out to be non-gradient. 
However, the homogenization tools facilitated the proof of the hydrodynamic limit, 
which made the proof much simpler than the usual proof of hydrodynamic limit for 
non-gradient processes (see for instance [5J Chapter 7]). 

We now describe the organization of the article. In Section [2] we define the W- 
Sobolev spaces and obtain some results, namely, approximation by smooth func- 
tions, Poincare's inequality, Rellich-Kondrachov theorem (compact embedding), 
and a characterization of the dual of the M^-Sobolev spaces. In Section[3]we define 
the VF-generalized elliptic equations, and what we call by weak solutions. We then 
obtain some energy estimates and use them together with Lax-Milgram's theorem 
to conclude results regarding existence, uniqueness and boundedness of such weak 
solutions. In Section[4]we define the VF-gcncralized parabolic equations, their weak 
solutions, and prove uniquenesses of these weak solutions. Moreover, a notion of 
energy is also introduced in this Section. Section [5j consists in obtaining discrete 
analogous results to the ones of the previous sections. This Section serves as pre- 
amble for the subsequent sections. In Section [6] we define the random operators 
we are interested and obtain homogenization results for them. Finally, Section [7J 
concludes the article with an application that is interesting for both probability 
and theoretical physics, which is the hydrodynamic limit for a gradient process in 
random environments. This application uses results from all the previous sections 
and provides a proof for existence of weak solutions of VT-gcneralized parabolic 
equations. 
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2. W-SOBOLEV SPACES 

This Section is devoted to the definition and derivation of properties of the 
VF-Sobolev spaces. We begin by introducing some notation, stating some known 
results, and giving a precise definition of these spaces in subsection 2.2. Subsection 
2.3 contains the proof of an approximation result. Poincare's inequality, Rellich- 
Kondrachov theorem and a characterization of the dual space of these Sobolev 
spaces are also obtained. 

Denote by T d = (M/Z) d = [0, l) d the d-dimcnsional torus, and by = 
(Z/NZ) d = {0, . . . , TV — l} d the d-dimensional discrete torus with N d points. 

Fix a function W : R d -> K such that 

d 

(3) W(x 1 ,...,x d )=J2w k (x k ), 

fc=i 

where each W k ■ K — * K is a strictly increasing right continuous function with left 
limits (cadlag), periodic in the sense that for all u G K 

W k (u + 1) - W k (u) = W k (l) - Wfe(0). 



Define the generalized derivative dw k of a function / : T d — > K by 

(4) 

o ,/ , ,. /(^lv ■ ;Xk + e, . . . ,x<j) - /(xi,. . .,X k ,. ■ .jXd) 
■ ;Xk, ...,x d ) = hm — — r -— r , 

e^o Wfc {xk+e)- W k (x k ) 

when the above limit exists and is finite. If for a function / : T d — > M. the generalized 
derivatives 9vK fc exist for all k, denote the generalized gradient of / by 

Vw/ = (^wi/, ■ ■ -,d Wd f) ■ 
Consider the operator Cw k ■ ^w k C L 2 (T) — > R given by 
(5) C W J = 9 Kfc d Wk /, 

whose domain is completely characterized in the following proposition: 
Proposition 2.1. The domain "Dw k consists of all functions f in L 2 (T) such that 

f(x) = a + bW k (x) + f W k (dy) [ j{z)dz 

J(0,x] JO 

for some function f in L 2 (T) that satisfies 

[ f(z)dz = and [ W k {dy)\b+ [ % ' f(z) dz) = 0. 
Jo J (o,i\ Jo 

The proof of Proposition 12.11 and further details can be found in [5] . Further- 
more, they also proved that these operators have a countable complete orthonormal 
system of eigenvectors, which we denote by A\y k - Then, following (IB] . 

d 

A w = {f -T l ^R:J(x u ... ] x d ) = l[f k (x k )J k eA Wk }, 



k=l 



where W is given by ©. 
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We may now build an operator analogous to £w k in T d . For a given set A, we 
denote by span(A) the linear subspace generated by A. Let H>w = span(Aw), and 
define the operator hw '■ h>w — * L 2 (T d ) as follows: for / = Ilfc=i fk G Aw, 

d d 

(6) h w {f){x 1 ,...x d ) = Y J fj 

k=lj=l,j^k 

and extend to 1D>vk by linearity. It is easy to see that if / € ©h/ 

d 

(7) L ff / = ^£wJ, 

fc=i 

where the application of C\y k on a function / : T d — > R is the natural one, i.e., it 
considers / only as a function of the fcth coordinate, and keeps all the remaining 
coordinates fixed. 

Let, for each k = 1, . . . , d, fk S Aw k be an eigenvector of Cw k associated to the 
eigenvalue Then / = Ilfc=i fk belongs to U>w and is an eigenvector of hw with 
eigenvalue £\=i Moreover, [TH] proved the following result: 

Lemma 2.2. The following statements hold: 

(a) The set B w is dense in L 2 (T d ); 

(b) The operator hw ■ hi>w — * L 2 (T d ) is symmetric and non-positive: 

(-h w f,f) > 0, 

where (•, •) is the standard inner product in L 2 (Y d ). 

2.1. The auxiliary space. Let L 2 k ^ Wk (T d ) be the Hilbert space of measurable 
functions H : T d — > R such that 

d{x k ®W k )H{x) 2 < oo, 

where d(x k ® Wk) represents the product measure in T d obtained from Lesbegue's 
measure in T d_1 and the measure induced by W k in T: 

d(x k ® Wk) = dx\ ■ ■ ■ dxk-i dWk dxk+i ■ • ■ dxd- 

Denote by (H, G) x k^ Wk the inner product of 1%^^ (T d ): 

(H,G) x ^ Wk = f d{x k ®W k )H{x)G{x) , 

and by || • || x fc^w - fc the norm induced by this inner product. 
Lemma 2.3. Let f,g€ ^>w> then for i = 1, . . . , d, 

(d Xl d w J(x))g(x) dx = - [ (d w J)(d Wz g)d(x l ^W l ). 



In particular, 



[ h w f(x)g(x) dx = -V / (d w J){d Wz g)d{x l ®Wi 

Jjd Jjd 
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Proof. Let /, g G Diy. By Fubini's theorem 



/ 



£wif{x)g{x)dx 



Jd-1 



C Wi f(x)g{x)dxi 



dx l 



where dx % is the Lebesgue product measure in T on the coordinates xi, . . . , 

Xj— 1, X^_(_i, . . . , Xrf. 

An application of [6l Lemma 3.1 (b)] and again Fubini's theorem concludes the 
proof of this Lemma. □ 



Let L\i^ w . (T d ) be the closed subspace of L 2 j(giW (T d ) consisting of the func- 
tions that have zero mean with respect to the measure d{x J ® Wj ) : 



fd{x j ®W 3 ) = 0. 



Finally, using the characterization of the functions in V>Wj given in Proposition 
12.11 and the definition of H>w, we have that the set {dw h\ h G ^w} is dense in 
(T d ). 



T 2 



2.2. The VF-Sobolev space. We define the Sobolev space of PF-generalized 
derivatives as the space of functions g G L 2 (T d ) such that for each i = l,...,d 
there exist fuctions Gi S L 2 i(g)W Q (T d ) satisfying the following integral by parts 
identity. 

(8) f {d Xi d w j)gdx = - f (d w J)G t d(x l ®W t ), 

for every function / G H>w- We denote this space by H\ t w^ d )- A stan- 
dard measure-theoretic argument allows one to prove that for each function g G 
H^ w {T d ) and i = 1, ...,d, we have a unique function Gi that satisfies (|8|). Note 
that B>w C H\ ) w^ d )- Moreover, if g G 3w then Gi = dwS- For this reason 
for each function g G H\ w we denote Gi simply by <9vV;Si an d we call it the ith 
generalized weak derivative of the function g with respect to W . 

Lemma 2.4. The set H ltW (T d ) is a Hilbert space with respect to the inner product 

d „ 

(9) (f,g)i,w = (f,9) + Y, / (dwJ)(d Wt g)d(x l ®W t ) 

Proof. Let (<7„) ra eN be a Cauchy sequence in Hi t w(T ), and denote by || • \\i,w the 
norm induced by the inner product |9]). By the definition of the norm || ■ Hi^w, 
we obtain that (g n )n£fi is a Cauchy sequence in L 2 (T d ) and that (9wi5n)neN is 
a Cauchy sequence in L 2 ^ w . (T d ) for each i = 1,.. .,d. Therefore, there exist 
functions g G L 2 (T d ) and Gi G L 2 x ^ w (T d ) such that g = lim^oo g n , and Gi = 
lim.,^00 dwign- It remains to be proved that Gi is, in fact, the ith generalized weak 
derivative of g with respect to W. But this follows from a simple calculation: for 
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each / £ U>w we have 



{d Xi dwif)gdx = lim / (d Xi d Wi f)g n dx 

lim / {d w J){d Wz g)d{x l ®Wi) 
(d w J)G t d(x l ®W t ), 

where we used Holder's inequality to pass the limit through the integral sign. □ 

2.3. Approximation by smooth functions and the energetic space. We will 
now obtain approximation of functions in the Sobolev space Hi t wC^ d ) by functions 
in H>w- Note that the functions in D^y can be seen as smooth, in the sense that 
one may apply the operator Lvy to these functions in the strong sense. 
Let us introduce (■, -}i_w the inner product on ~Bw defined by 

(10) (f,g)i,w = (f,g) + 

and note that by Lemma |2.3[ 

d 

i 



(f,9)i,w = {f, 9) + Y\ I (OwJKdwMx^Wi). 



Let Hi^wO?) De the set of all functions / in L 2 (T d ) for which there exists a 
sequence (/ n )neN m such that f n converges to / in L 2 (T d ) and /„ is a Cauchy 
sequence for the inner product (•, -)x,w- Such sequence (/ n )neN i s called admissible 
for /. 

For /, g in H 1>w (T d ), define 
(11) (f,g)i,W = lim (f n ,g n )i,w , 

n — >oc 

where (/„)„ e n, (g n )neN are admissible sequences for /, and g, respectively. By 
[T7| Proposition 5.3.3], this limit exists and does not depend on the admissible 
sequence chosen; the set l$w is dense in Hi y w] an d the embedding Hi y w C L 2 (T d ) 
is continuous. Moreover, HiwC^ d ) endowed with the inner product (v)i,w just 
defined is a Hilbert space. Denote H-j^ w the norm in H\ y w induced by (•, -)i,w- 
The space iJi ! v^(T d ) is called energetic space. For more details on the theory of 
energetic spaces see [TTl Chapter 5]. 

Note that H\ y w is the space of functions that can be approximated by functions 
in 3w with respect to the norm || • \\\,w- The following Proposition shows that this 
space is, in fact, the Sobolev space Hi t wC^ d )- 

Proposition 2.5 (Approximation by smooth functions). We have the equality of 
the sets 

In particular, we can approximate any function f in the Sobolev space iJx.wOT ) 
by functions in lD>w- 

Proof. Fix g G i/^vf (T d ). By definition, there exists a sequence g n in V>w such 
that g n converges to g in L 2 (T d ) and g n is Cauchy for the inner product (■, -)\,w- 
So, for each i = 1, . . . , d there exists functions Gi £ L^^Wi o0^ d ) sucn tliat dw^gn 
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converges to d in l? x ^ w . (T d ). Applying the Holder's inequality, we deduce that 
for every / G Byp 

(d Xi d Wi f)g dx = lim / (d Xl d Wl f)g n dx. 



By Lemma 12.31 wc obtain 

lim / (d Xi d w J)g n dx = lim / (d Wi f)(d Wt g n )d{x l ®W i ) 

n~ too Jj d n~ >oo Jj d 

(dwJ)Gi dix^W,). 

Then, g G Hi,w{T d ) and therefore H ltW (T d ) C #i, w (T d ). 

We will now prove that H\yy^ d ) is dense in Hi ! w(T d ), and since both of them 
are complete, they are equal. Note that since Oy^ is dense in L 2 (T d ) and Byy C 
#i,iv(T d ), we have that H l>w (T d ) is also dense in L 2 (T d ). 

Therefore, given a function g G _Hi,yt/-(T d ), w( 5 can approximate g by a sequence 
of functions {f n )nen in i/i,vv(T d ) with respect to the L 2 (T d ) norm. Let Fi >n be the 
ith generalized weak derivative of /„ with respect to W. We have, therefore, for 
each h G V>w 

lim / (d Wi h){F i<n -G i )d{x i ®W i ) = - lim / (d Xl d Wl h)(f n - g)dx = 0. 

n—>oo Jj, d ■ n->oo Jj d 

Denote by Ti, n : L xi ^ w . (T d ) — > R the sequence of bounded linear functionals 
induced by n — G^: 

for ft, G i^vK Q (T d ). We then note that, since the set {d\Yih;h G V>w} is dense 
in L 2 xi ^ w (T d ), Ti.n converges to pointwisely. By Banach-Steinhaus' Theorem, 
Ti^n converges strongly to 0, and, thus, F^ n converges to Gi in L xi ^ Wi (T d ), for 
each i = 1, . . . , d. Therefore, /„ converges to g in L 2 (T d ) and dwtfn converges to 
Gi in L 2 i<gjW (T d ) for each i, i.e., /„ converges to g with the norm || • \\i t w> an d 
the density of H hW (T d ) in H hW (T d ) follows. □ 

The next Corollary shows an analogous of the classic result for Sobolev spaces 
with dimension d = 1, which states that every function in the one-dimensional 
Sobolev space is absolutely continuous. 

Corollary 2.6. A function f in L 2 (T) belongs to the Sobolev space Hi : wC^) if and 
only if there exists F in L^-(T) and a finite constant c such that 

F(y)dW(y) = and f{x) = c + [ F(y)dW(y) 

(0,1] J(0,x] 

Lebesgue almost surely. 

Proof. In [6] the energetic extension H\yv(T) has the characterization given in 
Corollary 12.61 By Proposition 12.51 we have that these spaces coincide, and hence 
the proof follows. □ 
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From Proposition [231 we ma y use the notation Hiyy(T d ) for the Sobolev space 
Hi i wOF d )- Another interesting feature we have on this space, which is very useful 
in the study of elliptic equations, is the Poincare inequality: 



Corollary 2.7 (Poincare Inequality). For all f £ H\ t w(^ ) there exists a finite 
constant C such that 

2 

/ - / / dx 



„ 2 n ~ 

fdx < CV / (dwjfdix^Wi, 

JT d L 2 (T d ) 



C||V W /|| 2 



li^(T")- 

Proof. We begin by introducing some notations. For x,y £ T d , i = 0, . . . , d and 
teT, denote 

z(x,y,i) = (xi, . . . ...,y d )e T d 

and 

z(x, y, t, i) = (xi, . . . , x d -i, t, y d -i+2, ■ ■ ■ , yd) € T d . 
With this notation, we may write /(a;) — /(y) as the telescopic sum 

d 

f(x) - f{y) = ^2 f( z ( x ' V' i ~ !)) ~ 2/j *))• 



We are now in conditions to prove this Lemma. Let / £ 



then 



fdx\ 



< 



< 



T d 



j, i 



T d 



T d 



T d 



f(x) - f(y)dy 



L 2 (T d ) 

X / d w J(z(x,y,t,i))dW l (t)dy 

i=X J Vi 



d.r 



dx 



T d 



i=l 



1 2 



< c 



=i 

d 

E 

i=l 



E / I^VCi/CK^iM)) dWi{t)dy dx 

dwif{z(x,y,t,i)) dWd-i{t)® yd-i+i® ■■■®yd 

2 

d Wi f(z(x,y,t,i)) dW d -i{t)® dy d -t+i® ■ ■ -® dy d efe 



f/.C 



T d-i+l 



where in the next-to-last inequality, we used Jensen's inequality and the elementary 
inequality .t,) 2 < cc 2 for some positive constant C. To conclude the proof, 

one uses Proposition 12.51 to approximate functions in .Hi^CT^) by functions in 

□ 

2.4. A Rellich-Kondrachov theorem. In this subsection we prove an analogous 
of the Rellich-Kondrachov theorem for the VF-Sobolev spaces. We begin by stating 
this result in dimension 1, whose proof can be found in [6l Lemma 3.3]. 

Lemma 2.8. Fix some k £ {l,...,d}. The embedding Hiw k (T) C L 2 (T) is 
compact. 



10 



A.B. SIMAS AND F.J. VALENTIM 



Recall that they proved this result for the energetic extension, but in view of 
Proposition [231 this result holds for our Sobolev space H\ y w k i^)- 

Proposition 2.9 (Rellich-Kondrachov). The embedding H lM {T d ) C L 2 (T d ) is 
compact. 

Proof. We will outline the strategy of the proof. Using the definition of the set H>w 
and the fact that it is dense in Hi t w{^ d ), it is enough to show this fact for sequences 
in Dvy. From this point, the main tool is Lemma [2751 and Cantor's diagonal method 
to obtain converging subsequences. 

We begin by noting that by Proposition ^. 5[ it is enough to prove that the embed 
D w C L 2 (T d ) is compact. 

Let C > and consider a sequence {v n ) n ^jq in D^/, with 1 1 w T1 1 j i , < C for all 
n £ N. We have, by definition of V>w (see the definition at the beginning of Section 
H}, that each v n can be expressed as a finite linear combination of elements in Aw- 
Furthermore, each clement in Aw is a product of elements in Aw k for k — 1, . . . , d. 
Therefore, we can write v n as 

N(n) d N(n) 

3=1 k=l 3=1 

where g^j £ Aw k ,oij 6R, = IIfc=i Sj'k' an d N(n) is chosen such that N(n) > n 
(we can complete with zeros if necessary). Recall that these functions g% ■ have 
HffrjlU 2 (T) = 1, and hence, ||gJ||z,2 (T «£ ) = 1. Moreover, the set {g?, . . . , (n) } is 
orthogonal in L 2 (T d ). 

From orthogonality, we obtain that 

N(n) 

E «) 2 < C 2 , uniformly in n E N. 

3 = 1 

Note that the uniform boundedness of v n in Hi_w(^ d ) implies the uniform 
boundedness of ||.9fej||i,w fc j for all k = 1, . . . , d, j = 1, . . . , N(n) and n € N. Our 
goal now is to apply Lemma 1^51 to our current setup. 

Consider the sequence of functions 1 in C^) ■ By Lemma [2.81 this 

sequence has a converging subsequence, and we call the limit point ai<7i,i. Repeat 
this step d— 1 times for the sequences g - ^ 1 in H\.w h (T), for k = 2, . . . , d, considering 
in each step a subsequence of the previous step, to obtain converging subsequences, 
and call their limit points gk,i- At the end of this procedure, we obtain a converging 
subsequence of rifc=i a i9ik> w ^th limit point rife=i a i5i,fe € L 2 (T d ), which we will 
denote by u\g\. 

In the jth step, in which we want to obtain the limit point ct^g^ we repeat the 
previous idea, with the sequences a^g^ and g™ k , with n < j and k = 2, . . . , d. We 
note that it is always necessary to consider a subsequence of all the previous steps. 

This procedure provides limiting functions otjgj, for all j G N. From now on, 
we use the notation v n to mean the diagonal sequence obtained to ensure the 
convergence of the functions o^g™ to otjgj. We claim that the function 

oo 

v = E "j'l. 
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is well-defined and belongs to L 2 (T d ). To prove this claim, note that the set {gfcjfceN 
is orthonormal by the continuity of the inner product. Suppose that there exists 
N G N such that 

f>,) 2 > c 2 - 

3 = 1 

We have that the sequence of functions 

N 
3=1 

converges to 

N 

N \ " 

3 = 1 

Since ||v^|| < C uniformly in n G N, this yields a contradiction. Therefore v G 
L 2 {T d ) with the bound ||u|| < C. 

It remains to be proved that v n has a subsequence that converges to v. Choose 
N so large that ||w -v N \\ < e/3, -v N \\ < e/3 and \\v% - v„|| < e/3, and use 
the triangle inequality to conclude the proof. 

□ 

2.5. The space H^(T d ). Let H^(T d ) be the dual space to H hW (T d ), that is, 
if^r 1 (T d ) is the set of bounded linear functional on Hi^ w (T d ). Our objective in 
this subsection is to characterize the elements of this space. This proof is based on 
the characterization of the dual of the standard Sobolev space in R d (see [2]). 
We will write (•, •) to denote the pairing between H^ l (T d ) and H hW (T d ). 

Lemma 2.10. / G H^(T d ) if and only if there exist functions fo G L 2 (T d ), and 
fk e Ll 

d 

(12) f = fo-J2 d *J" 

i=l 

in the sense that for v G -Hi,w(T d ) 

(/,«)=/ hvdx + Y^f Md Wl v)d(x l ®W t ). 

Furthermore, 




Proof. Let / G H w (T d ). Applying the Riesz Representation Theorem, we deduce 
the existence of a unique function u G i/i^CI" 11 ) satisfying (/, v) = (u, v)i i w, for 
all v G i?i^(T d ), that is 

(13) f uvdx + Y j [ {d W] u){d W] v)d{x3®W 3 ) = (/,«), for all « e Fi, w (T d ). 



12 



A.B. SIMAS AND F.J. VALENTIM 



This establishes the first claim of the Lemma for fo = u and fi = dwiU, for 
i = 1 , . . . , d. 

Assume now that / £ H^, (T"), 



(14) 



JT d i=l ^ Td 



for go, 5i, ■ ■ ■ ,9d € L 2 x](g)Wj (T d ). Setting v = u in |T3]), using JTU), and applying 
the Cauchy-Schwartz inequality twice, we deduce 

d 

(15) h 2 



\i,w< 9odx + J2 dw.gfdix^Wi). 



From (|T3")) it follows that 
if HfH i,w < 1- Consequently 

II/IIh- < ||«||i lW r. 

Setting u = in ()13|) . we deduce that, in fact, 

The result now follows from the above expression and equation (|15|) . □ 

3. W- Generalized elliptic equations 

This subsection investigates the solvability of uniformly elliptic generalized par- 
tial differential equations defined below. Energy methods within Sobolev spaces 
are, essentially, the techniques exploited. 

Let A = (au(x))dxd, x £ T d , be a diagonal matrix function such that there exists 
a constant 6 > satisfying 

(16) 6- 1 <a u (x)<6, 

for every x £ T d and i = 1, . . . , d. To keep notation simple, we write a, (a;) to mean 
a u (x). 

Our interest lies on the study of the problem 

(17) T\u = /, 

where u : T d — > M. is the unknown function and / : T d — > K is given. Here T\ 
denotes the generalized elliptic operator 

d 

(18) T\u := Aw — Vj4Vw-u := Ait — 5^ (ai(x)dwiu) ■ 

i=l 

The bilinear form B[-, ■] associated with the elliptic operator T\ is given by 

d . 

(19) B[u,v] = X(u,v) a.i{x)(d Wt u)(d Wz v) d{Wi®x t ), 

i=l •* 

where u, v £ Hi.yy(T d ). 

Let / £ H^-(T d ). A funct ion u £ Hi ! w{'K d ) is sa id to be a weak solution of the 
equation T\u = / if 

B[u,v] = (f,v) for all v £ ff w (T d ). 
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Recall a classic result from linear functional analysis, which provides in certain 
circumstances the existence and uniqueness of weak solutions of our problem, and 
whose proof can be found, for instance, in [2]. Let TL be a Hilbert space endowed 
with inner product < ■, ■ > and norm ||| • |||. Also, (•, •) denotes the pairing of TL 
with its dual space. 

Theorem 3.1 (Lax-Milgram Theorem). Assume that B : TL x TL — > R is a bilinear 
mapping on Hilbert space TL, for which there exist constants a > and (3 > such 
that for all u 1 v E TL, 

\B[u,v]\ < a\\\u\\\ ■ \\\v\\\ and B[u,u] > f3\\\u\\\ 2 . 

Let f : TL — *■ R be a bounded linear functional on TL. Then there exists a unique 
element u £ TL such that 

B[u,v] = (f,v), 

for all v E TL. 

Return now to the specific bilinear form B[-,-} defined in (fH)|) . Our goal now is 
to verify the hypothesis of Lax-Milgram Theorem for our setup. We consider the 
cases A = and A > separately. We begin by analyzing the case in which A = 0. 

Let H^- w (Y d ) be the set of functions in Hi ^(T d ) which arc orthogonal to the 
constant functions: 

Htw(^ d ) = {/ 6 / fdx = 0}. 

The space H^- w (T d ) is the natural environment to treat elliptic operators with 
Neumann condition. 

Proposition 3.2 (Energy estimates for A = 0). Let B be the bilinear form on 
Hi yy(T d ) defined in (|19p with A = 0. There exist constants a > and j3 > such 
that for all u, v E H\^wi^ d ), 

\B[u,v}\ < a\\u\\ ltW \\v\\i t w 

and for all u E H^ w 

B[u,u] >P\\u\\l w . 

Proof. By (|16[) . the computation of the upper bound a easily follows. For the lower 
bound /3, we have for u E H^ w (T d ), 

\\u\\lw= I u 2 dx + Y^I (d Wl u) 2 d(x l ®Wi). 

Using Poincare's inequality and (fT6|) . we obtain a constant C > such that the 
previous expression is bounded above by 

C j (d Wz u) 2 d(x l ® W t ) < CB[u, u]. 

The lemma follows from the previous estimates. □ 

Corollary 3.3. Let f E L 2 (T d ). There exists a weak solution u E Hi. w (T d ) for 
the equation 

(20) VAV w u = f 
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if and only if 

[ fdx = 0. 

In this case, we have uniquenesses of the weak solutions if we disregard addition by 
constant functions. Also, let u be the unique weak solution of (|20p in H^ w (T d ). 
Then 

||w||i,W < CII/IIl 2 ^), 
for some constant C independent of f . 

Proof. Suppose that there exists a weak solution u £ H\^w^ d ) of (|20p . Since the 
function v = 1 £ Hi t wC^ )> we have by definition of weak solution that 

/ fdx = B[u,v] = 0. 

Now, let / £ L 2 (T d ) with J Jd fdx = 0. Consider the bilinear form B, defined in 
(Unj) with A = 0, on the Hilbert space H^ w (T d ). By Proposition [XI B satisfies 
the hypothesis of the Lax-Milgram's Theorem. Further. / defines the bounded 
linear functional in H^ w (T d ) given by (/,<?) = (/, g) for every g £ H^ w (T d ). 
Then, an application of Lax-Milgram's Theorem yields that there exists a unique 
u £ H^ w {T d ) such that 

B[u,v] = (f,v) for all v £ #iV( T<i )- 
Moreover, by Proposition 13.21 there is a (3 > such that 

- = (/)") ^ ll/IU 2 (Ti)||'"||i=(T'2) < ll/lli2(T'»)||«||l,W- 

The existence of weak solutions and the bound C in the statement of the Corol- 
lary follows from the previous expression. □ 

We now analyze the case in which A > 0. 

Proposition 3.4 (Energy estimates for A > 0). Let f £ L 2 (T d ). There exists a 
unique weak solution u £ £fi i w(T d ) for the equation 

(21) \u-S7AV w u = /, A > 0. 

This solution enjoys the following bounds 

\\u\\i,w < C||/|| L 2 (T d ) 

for some constant C > independent of f , and 

\\U\\ < A _1 ||/||£2(Ti). 

Proof. Let (3 = min{X, O^ 1 } > and a = max{\, 9} < oo, where 8 is given in (|16|) . 
An elementary computation shows that 

< a||tt||i,w \\v\\i,w and B[u,u] > /9||w||i iW -- 

By Lax-Milgram's Theorem, there exists a unique solution u £ i/i^ (T d ) of ((2"Tj) . 
Note that 

^ = (/,U) < ||/|| i 2 (T d ) ||u|| L 2 (T£i) < ||/|| i 2 (T d ) ||u||i :V i/, 

and therefore < C||/ll,L 2 (T d ) f° r some constant C > independent of /. 

The computation to obtain the other bound is analogous. □ 
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Remark 3.5. Let : B w -> L 2 (T d ) be given by Ljy = VAVw . TTws operator 
has the properties stated in Theorem 2. 1 in [16j . We now outline the main steps to 
prove it. Following |16| . we may prove an analogous of Lemma \2.2\ for the operator 
L^. Using the bounds on the diagonal matrix A and Provosition \2.9\ (Rellich- 
Kondrachov), we conclude that the energetic extension of the space induced by this 
operator has compact embedding in L 2 (T d ). The previous results together with \17\ 
Theorems 5. 5. a and 5.5.c] implies that has a self-adjoint extension Cyy, which 
is dissipative and non-positive, and its eigenvectors form a complete orthonormal 
set in L 2 (T d ). Furthermore, the set of eigenvalues of this extension is countable 
and its elements can be ordered resulting in a non-increasing sequence that tends to 
— oo. 

Remark 3.6. Let £^r be the self-adjoint extension given in Remark \3.5l and T>^ 
its domain. For A > the operator AI — Cyy : T>w ~ * L 2 (T d ) is bijective. Therefore, 
the equation 

Xu - VA\7 w u = f, 

has strong solution in D^y if and only if f G (AI — £^,,)(Dw), where I is the 
identity operator and (AI— £^)(Ojy) stands for the range o/Bjy under the operator 
AI — C^r . Moreover, this strong solution coincides with the weak solution obtained 
in Proposition \3.4\ 

4. W- Generalized parabolic equations 

In this Section, we study a class of VT-gcneralizcd PDEs that involves time: the 
parabolic equations. The parabolic equations are often used to describe in physical 
applications the time-evolution of the density of some quantity, say a chemical 
concentration within a region. The motivation of this generalization is to enlarge 
the possibility of such applications, for instance, these equations may be used to 
model a diffusion of particles within a region with membranes (see [6l 116]). 

We begin by introducing the class of VF-generalized parabolic equations we are 
interested. Then, we define what is meant by weak solution of such equations, using 
the VK-Sobolev spaces, and prove uniquenesses of these weak solutions. In Section 
we obtain existence of weak solutions of these equations. 

Fix T > and let (B, || ■ be a Banach space. Wc denote by L 2 ([0, T],E>) the 
Banach space of measurable functions U : [0, T] — ► B for which 

\\U\\h {[0 .T],B) ■■= f T \\U t \\ 2 B dt < ^. 

Jo 

Let A = A(t, x) be a diagonal matrix satisfying the ellipticity condition (fT6|) for 
all t € [0,T], $ : [l,r] — > R be a continuously differentiable function such that 

B- 1 < $'(2;) < B, 

for all x, where B > 0, l,r € K are constants. We will consider the equation 

f d t u = VAV w ^(u) in (0,T] x T d , 
{ > \ u = 7 in {0} x T d . 

where u : [0, T] xr J -tl is the unknown function and 7 : T d — > R is given. 
We say that a function p = p(t, x) is a weak solution of the problem ([22"]) if: 
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• For every H E V>w the following integral identity holds 
p(t,x)H(x)dx - [ j(x)H(x)dx = [ [ <&{p{s , x))V AV w H (x)dx ds 

JT d JO Jf d 

• $(p(-, •)) and p{; •) belong to L 2 ([0, T],H 1 , w (T d )): 

pT 

Mp{s,x))\\l HTd) + \\V w $(p(s,x))\\l 2 {Td) ds < 00, 

./ u 

and 



Consider the energy in jth direction of a function u(s, x) as 

Q./(lt) = sup <2 / (d Xj dw j H)(s,x)u(s,x)dxds 
H£Dw L Jo JT d 

- f ds f [dw^is^tfdix^Wj)}, 
Jo Jr d ' 

and the total energy of a function u(s, x) as 

d 

The notion of energy is important in probability theory and is often used in large 
deviations of Markov processes. We also use this notion to prove the hydrodynamic 
limit in Section[7J The following lemma shows the connection between the functions 
of finite energy and functions in the Sobolev space. 

Lemma 4.1. A function u E L 2 ([Q,T], L 2 (T d )) has finite energy if and only if u 
belongs to L 2 ([0,T], Hi } wO? d ))- In the case the energy is finite, we have 

Q(u) = J^ \\Vwu\\ 2 L 2 w{Td) dt. 

Proof. Consider functions U E L 2 ([0,T], L 2 j( ^ w Q (T d )) as trajectories in 
Ll j(giW . (T d ), that is, consider a trajectory U : [0,7] -> L 2 x j^w 3 .o^ d ) and denne 
U(s,x)as U{s,x) := [U(s)](x). 

Let u E L 2 ([0,T],L 2 (T d )) and recall that the set {d W] H:H E B w } is 
dense in L 2 x ^ w , Q (T d ). Then the set {d Wj H{s,x); H G L 2 ([0,T],D w )} is dense 
in L 2 ([0,T], L 2 j ^ )W (T d )). Suppose that u has finite energy, and let H E 
L 2 {[0,T],D w ), then 

Tj{d W] H) = I I (d Xj d Wj H){s,x)u(s,x)dxds 
Jo Jr d 

is a bounded linear functional in L 2 ([0, T] , L 2 jgiW . (T d )). Consequently, by Riesz's 
representation theorem, there exists a function Gj E L 2 ([0,T], L 2 ^ w (T d )) such 
that 

Tj(dw 3 H) — I I (dw j H)(x)G j (s,x)dxds, 
Jo Jr d 

for all H E L 2 {[0,T],D w ). 
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From the uniqueness of the generalized weak derivative, we have that Gj(s, x) = 
-dWjU{s,x). 

Now, suppose u belongs to L 2 ([0, T], H liW (T d )) and let H £ L 2 ([0, T],B W ). 
Then, we have 

2/ / (d Xj d W] H)(s,x)u(s 1 x)dxds- [ dsf (d Wj H(s,x)) 2 d(x j ® Wj) = 

JO JT d JO JT d 

-2 f [ d W -H(s,x)d Wj u(s,x)d{x j ig)W j )- f [ (d Wj H(s,x)) 2 d(x j ®Wj) 
Jo Jt<* Jo Jt<* 

We can rewrite the right-hand side of the above expression as 

(23) - 2(d Wj H, 2d Wj u + d Wj H) xmw . . 

A simple calculation shows that, for a Hilbcrt space TL with inner product <■, ■>, 
the following inequality holds: 

1 

— <v. u + v> < — <u, u>, 
~ 4 

for all u, v £ 7i, and we have equality only when v = —l/2u. 
Therefore, by the previous estimates and (|23l) 

(d Xj dw j H)(s,x)u(s,x)dxds — / ds / [pWj H{s, x)) 2 d(x J : ® Wj) < 
Ji d Jo Jt d 

T r 

I (d W] u{s,x)) 2 d(x j ® Wj). 

I0 JT d 

By the definition of energy we have for each j = 1, . . . , d, 

Qj(u)< [ I (dw.u(s,x)) 2 d{x j ®Wj). 
Jo Jr d 

Hence, the total energy is finite. Using the fact that £ 2 ([0, T], Dvy) is dense in 
L 2 ([0,T],H hW (T d )), we have that 



A 1 " 

IIVm/uH 2 ^ (jd^dt. 



□ 

4.1. Uniqueness of weak solutions of the parabolic equation. Recall that we 
denote by (•, •) the inner product of the Hilbert space L 2 (T l ). Fix H,G £ L 2 (T d ), 
A > 0, and denote by H\ and G\ in H^wC^ ) the unique weak solutions of the 
elliptic equations 

\H\ - VAV W H X = H, 

and 

XG X - VA\7 W G X = G, 
respectively Then, we have the following symmetry property 

(G X ,H) = (G,H X ). 
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In fact, both terms in the previous equality are equal to 

A / H x Gx+y2a n [ (d Wj H x )(d Wj Gx)d(xi ®Wj). 



3 = 1 



Let p : K+ xT^ [l,r] be a weak solution of the parabolic equation (J^U). Since 
p, $0) 6 L 2 ([0,T],H 1>w (T d )), and the set B w is dense in H 1>w (T d ), wc have for 
every H in H hW {T d ), ' 



(24) ( Pt ,H) - < 7l tf) 



y2" !! / { d w j ^{p s ),dw s H) x i 9Wj ds 

3 = 1 Jo 



for all t > 0. 

Denote by p x £ Hi } wC^ d ) the unique weak solution of the elliptic equation 

(25) Xp x s - VAV w p x = P (s, ■). 
We claim that 

(26) (pt ,Pt) - (po,Po) = -2^ a « / (dwj^iPs) , dw 3 Ps)x^w 3 ds 

3=1 J ° 

for all t > 0. 

To prove this claim, fix t > and consider a partition = to < ti < • • • < t n = t 
of the interval [0,t]. Using the telescopic sum, we obtain 

n-l 

(Pt , Pt) - (Po , Pq) = J2(Pt k+1 , Pt k+1 ) - {pt k+1 , Pt k ) 

n-l 

+ ^2(pt k+1 , Pt k ) - <pt fc . Pt k ) ■ 

k=Q 

We handle the first term, the second one being similar. From the symmetric 
property of the weak solutions, p x k 1 belongs to £?i ; w(T d ) and since p is a weak 
solution of (1221), 



d rtk+i 

{Pt k +i , Pt k+1 ) - (Pt k+1 ,Pt k ) = ~^Z a 33 I (dw^iPs), d Wj Pt k+1 ) ds ■ 

3 = 1 Jtk 

Add and subtract (9w.$(p s ), dw/jP X ) inside the integral on the right hand side 
of the above expression. The time integral of this term is exactly the expression 
announced in (|26p and the remainder is given by 

J2 a n I {(dwMPs) > dwjPs) ~ (d w ^( Ps ) , d W]P x k+1 )) ds . 

3 = 1 tk 

Since p x is the unique weak solution of the elliptic equation (f2"5")) . and the weak 
solution has the symmetric property, we may rewrite the previous difference as 

- (Hp s ),p s )} - \{(<f( P s)\ Ptk+1 ) - mp s )\ps)} . 

The time integral between and tk+i of the second term is equal to 
-A / " +1 ds f (d Wj ^( Ps ) x , dw^(p r )) dr 
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because p is a weak solution of (|22|) and &(p s ) belongs to H\ y wi^ d ). It follows 
from the boundedness of the weak solution given in Proposition 13.41 and from the 
boundedness of the L 2 je>w (T d ) norm of dwj^(p) obtained in expression (f2"4"|) . that 
this expression is of order (tk+i — tk) 2 . 

To conclude the proof of claim (|26|) it remains to be shown that 

?±Z} fth+i , , 

£ / {{Hp s ),Pt k+1 ) - m Ps ), Ps )\ds 

k=0 Jtk 

vanishes as the mesh of the partition tends to 0. Using, again, the fact that p is a 
weak solution, we may rewrite the sum as 

-J2 f k+1 ds f ^ (d Wj <S>(p s ) , dw.$(p r )) dr . 

k=0 Jtk Js 

We have that this expression vanishes as the mesh of the partition tends to 
from the boundedness of the L 2 ^ w . (T d ) norm of 3^ This proves (f2l)|) . 

Recall the definition of the constant B given at the beginning of this Section. 

Lemma 4.2. Fix X > 0, two density profiles j , 7 2 : T — > [l,r] and denote by p , 
p 2 weak solutions of (|22[1 with initial value 7 , 7 2 , respectively. Then, 



(pl-p 2 ,pl<^p^) < ( 7 i- 7 2 )7 M_ 7 ^ e ^/ 2 
for all t > 0. In particular, there exists at most one weak solution of 
Proof. We begin by showing that if there exists A > such that 

(H, H x ) — 0, 
then H = 0. In fact, we would have the following 

f A(i? A ) 2 dx + Vajj [ (d Wj H x f d{x j ®W 3 )= f HH x dx = 0, 

JT d j =1 JT d JT d 

which implies that ||i? A || ffl „,(T d ) = 0, and hence H\ = 0, which yields H = 0. 

Fix two density profiles 7 1 , 7 2 : T d — > [l,r). Let p 1 , p 2 be two weak solutions 
with initial values 7 1 , 7 2 , respectively. By (|26|) . for any A > 0, 



h^{p\)-<5>{p 2 s ),pl-p 2 s )ds + 2Xj^(<P(pl)-<P(p 2 ),p^-p 2 s x )ds. 

Define the inner product in Hi : w{^ d ) 

(u,v) x = (u,v x ). 

This is, in fact, an inner product, since (u, v)\ = {v, u)\ by the symmetric property, 
and if u 7^ 0, then {u, u}\ > 0: 



uu\dx = X 



/ u\dx + Y] ajj / (d Wj u x ) 2 d(x 3 ® Wj). 

JT d ~~f JT d 



The linearity of this inner product can be easily verified. 
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Then, we have 

2\J* - Hp 2 s ) , P l< X pf X ) ds = 2X^ - $(p 2 s ) , p] - fi)ds. 

By using the Cauchy-Schwartz inequality twice, the term on the right hand side 
of the above formula is bounded above by 

1 J* (HpD , $(pir - $(P 2 S ) X ) ds + AX 2 J* (pi pi , fi* - ds 

for every A > 0. From Proposition 13. 4i we have that ||w A || < A _1 ||u||, and since 
is bounded by B, the first term of the previous expression is less than or equal to 

■^J* (pl-p 2 s ,S>(pl)-<S>(p*))ds. 

Choosing A = B/2X, this expression cancels with the first term on the right hand 
side of (|27|) . In particular, the left hand side of this formula is bounded by 

^ l (pl - Pi Pl- X - pf X ) ds . 

To conclude, recall Gronwall's inequality. □ 

Remark 4.3. Let C^/ : T>w — > L 2 (T d ) be the self-adjoint extension given in Re- 
mark [375[ For X > 0. define the resolvent operator Gf = (XI — Cy^) 1 - Following 
[6),I16|. another possible definition of weak solution of equation (|22|) is given as fol- 
lows: a bounded function p : [0, T] x T d [l,r] is said to be a weak solution of the 
parabolic differential equation (|22p if 

(28) (fH,G£h) - (l,Gih) = [ (§{p s ),L^Gih)ds 

Jo 

for every continuous function h : T d — > R, t £ [0, T], and all X > 0. We claim 
that this definition of weak solution coincides with our definition introduced at the 
beginning of Section^ Indeed, for continuous h : T d — > R, G^h belongs to T>w- 
Since Bty is dense in T>w with respect to the Hi\y(T d )-norm, it follows that our 
definition implies the current definition. Conversely, since the set of continuous 
functions is dense in L 2 (T d ), the identity (|28p is valid for all h £ L 2 (T d ). Therefore, 
for each H £ T>w we have 

(p t ,H) - (j,H) = f (<S>(p s ),£*H)ds. 
Jo 

In particular, the above identity holds for every H £ and therefore the 

integral identity in our definition of weak solutions holds. 

It remains to be checked that the weak solution of the current definition belongs 
to L 2 ([Q, T], Hi^w(T )). This follows from the fact that there exists at most one 
weak solution satisfying (|28[) . that this unique solution has finite energy, and from 
Lemma \4-l\ A proof of the fact that there exists at most one solution satisfying 
(|28p, and that this unique solution has finite energy, can be found in [6lll6j. 

Finally, the integral identity of our definition of weak solution has an advantage 
regarding the integral identity (|28|) , due to the fact that we do not need the resolvent 
operator G^ for any X. Moreover, we have an explicit characterization of our test 
functions. 
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5. VF-Generalized Sobolev spaces: Discrete version 



We will now establish some of the results obtained in the above sections to the 
discrete version of the H^-Sobolev space. Our motivation to obtain these results is 
that they will be useful when studying homogenization in Section [6l We begin by 
introducing some definitions and notations. 

Fix W as in ((3]) and functions /, g defined on A r_1 T^ r . Consider the following 
difference operators: d^f. , which is the standard difference operator, 



N 



f 



N 



f 



and dy!r. , which is the Wj-difference operator: 



/(-^) /(*) 



for x € Tfj. We introduce the following scalar product 



(f,9)N ■= 

x£T d N 

^-,N :=^rr E mg(x)(W((x + ej )/N)-W(x/N)), 



(f,g)i,w,N ■=(/, 9)n + ^2(dwJ,d^g) w ^ N , 

3 = 1 

and its induced norms 



- (/i/)jv, II /II i,^ (f^.) - (f>f)Wj,N and \\f\\ Hlw (jd) - (f,f)i,w,N- 



These norms are natural discretizations of the norms introduced in the previous 
sections. Note that the properties of the Lebesgue's measure used in the proof of 
Corollary 12.71 also holds for the normalized counting measure. Therefore, we may 
use the same arguments of this Corollary to prove its discrete version. 

Lemma 5.1 (Discrete Poincarc Inequality). There exists a finite constant C such 
that 



f 



N d 



x£T d 



< c||v^/|u (Tf{ 



L 2 (Tf 



whe 



for all f : N^T 



\ V wf\\l lATlv) 



£lK-/l 



-lTd 

N 



Let A be a diagonal matrix satisfying |T 
problem 



We are interested in studying the 



(29) 
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where u : N — > R is the unknown function, / : N — > R is given, and 



denotes the discrete generalized elliptic operator 

rj*u := Xu-V N A V] 



(30) rfu := \u-W N AV^u, 



with 

V W AV^ U := J2dZ( ai (x/N)d^u 

i=l 

The bilinear form B N [, ■] associated with the elliptic operator Tj^ is given by 

B N [u,v] = \{u,v) N + 

(31) 

where u,v : iV _1 T^ -> R. 

A function u : N^ l T c ^ — » R is said to be a weak solution of the equation 
T»u = f if 

?-?Vr„. „,i _ / -f „,\__ „n „, . a;-- l-ipd 



B JV [u,w] = (f,v) N for all u : N~ T% 



Wc say that a function / : iV 1 Tf f — > R belongs to the discrete space of functions 



orthogonal to the constant functions H^(Tff) if 

± y: hx/n) = o. 

The following results are analogous to the weak solutions of generalized elliptic 
equations for this discrete version. We remark that the proofs of these lemmas are 
identical to the ones in the continuous case. Furthermore, the weak solution for the 
case A = is unique in H^(Tff). 

Lemma 5.2. The equation 

V N AV^u = /, 

has weak solution u : i\ r_1 T^ r — > R if and only if 



w E /(*) = °- 



N d 

In this case we have uniqueness of the solution disregarding addition by constants. 
Moreover, if u € Hj^(T%) we have the bound 

\H\hi, w {t%) < C \\f\\i?<j%)> and < ^^l/IL^x^), 
where C > does not depend on f nor N . 

Lemma 5.3. Let A > 0. There exists a unique weak solution u : 7V~ 1 T'^r — > R of 
the equation 

(32) \u - W N AV^u = f. 

Moreover, 

hllXiMK) ~ C Wf\\LZ(T%), and ll«lli=>(T*) < ^ll/lll^T*), 

where C > does not depend neither on f nor N . 
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Remark 5.4. Note that in the set of functions in Tfj we have a "Dirac measure" 
concentrated in a point x as a function: the function that takes value N d in x and 
zero elsewhere. Therefore, we may integrate these weak solutions with respect to 
this function to obtain that every weak solution is, in fact, a strong solution. 

5.1. Connections between the discrete and continuous Sobolev spaces. 

Given a function / G iJi j vi/(T' i ), we can define its restriction /jv to the lattice 
N- l T d N as 

f N {x) = f(x) if x e N^T d N . 

However, given a function / : N~ 1 T d v — ► R it is not straightforward how to 
define an extension belonging to To do so, we need the definition of 

M^-interpolation, which we give below. 

Let $n ■ N'^-Tn — > R and W : R — > M, a strictly increasing right continuous 
function with left limits (cadlag), and periodic. The I^-intcrpolation f^ of /jv is 
given by: 

r( , W«x + l)/N)-W((x + t)/N) 

t^ X + t > - W((x + 1)/N)-W(x/N) }{X) + 



W((x + t)/N)-W(x/N) 

W((x + 1)/N) - W{x/N) J{ ' 



for < t < 1. Note that 



9fk (x+t) _ f(x + l)-f(x) _ N 

dW { + >~ W((x + 1)/N)-W(x/N) ~ ° wJ[x> - 

Using the standard construction of d-dimensional linear interpolation, it is pos- 
sible to define the VF-interpolation of a function /jv : — > R, with W(x) = 
J2i=i Wi(xi) as defined in (|3|). 

We now establish the connection between the discrete and continuous Sobolev 
spaces by showing how a sequence of functions defined in Tfj can converge to a 
function in H\^w^ d )- 

We say that a family f^ G L 2 (T^) converges strongly (rcsp. weakly) to the 
function / G L 2 (T d ) as N — ► oo if f^ converges strongly (resp. weakly) to the 
function /. From now on we will omit the symbol " * " in the VF-interpolated 
function, and denoting them simply by f^. 

The convergence in H^(T d ) can be defined in terms of duality. Namely, we say 
that a functional /at on Tfj converges to / G H^~(T d ) strongly (resp. weakly) if 
for any sequence of functions un : Tfj — * R and u G Hi t wC^ d ) such that un — > u 
weakly (resp. strongly) in if^wCF*), we have 

Un,u n ) n — >{f,u), as N -> oo. 

Remark 5.5. Suppose in Lemma \5.3\ that f G L 2 (T d ), and iet u be a weak solution 
of the problem (|32[) , £/ien we /iai>e i/ie following bound 

II u IIj?i,w(t^) - C'll/lli 2 ^): 
since H/H^t* ) -> ||/|| L 2(T«) as iV ^ oo. 
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6. HOMOGENIZATION 



In this "brief" Section we prove a homogenization result for the W^-gcneralizcd 
differential operator. We follow the approach considered in [TI]. The study of ho- 
mogenization is motivated by several applications in mechanics, physics, chemistry 
and engineering. The focus of our approach is to study the asymptotic behavior 
of effective coefficients for a family of random difference schemes whose coefficients 
can be obtained by the discretization of random high-contrast lattice structures. 

This Section is structured as follows: in subsection 6.1 we define the concept of H- 
convergence together with some properties; subsection 6.2 deals with a description 
of the random environment along with some definitions, whereas the main result is 
proved in subsection 6.3. 

6.1. ii-convergence. We say that the diagonal matrix A N = (ajj) //-converges 

to the diagonal matrix A = {ajj), denoted by A N A, if, for every sequence 
f N £ H^(T d N ) such that /-^ ->• / as iV ->• oo in H^(T d ), we have 

• un — > uo weakly in H\,wi^ d ) as N — > oo, 

• afjd^UN -> ajjdwjUo weakly in L\ j%w . (T d ) for each j = 1, . . . , d, 
where un ■ Tj^ — > M is the solution of the problem 



The notion of convergence used in both items above was defined in subsection 15. II 
We now obtain a property regarding if- convergence. 



where f G H w L (T d ) is fixed. Then, the following limit relations hold true: 



and uo S Hi^ w (T d 



\u N -V N A N S7^u N = f N 
) is the solution of the problem 
Amo — VAVwUq = f . 



Proposition 6.1. Let A 





and 



1 



d 



N d-i 



EE 



a^{x){d^u N {x)) 2 [WjUxj + 1)/N) - Wjixj/N)} 




as N 



oo. 



Proof. We begin by noting that 



(33) 
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as N — > oo since un — uq converges weakly to in iJi,w(T d ). On the other hand, 
we have 

j^d E f( UN ~ u °) = Y (^un -V N A N V^u N )(u N -u Q ) 

— iv . i\r_ .jv 



E uj\[Uq + —J E u \7 N A N \7? w z 



N d ^ v " u 1 jv d 

Using the weak convergences of mtv and ajjdyy.UN, and the convergence in (|33p . we 
obtain, after a summation by parts in the above expressions, 



(34) N -^?\l utdx + V / a^^uo)*^®^). 



JT d - =l JT d 



By Lemma |5.3[ the sequence un is || • \\i.w bounded uniformly. Suppose, now, that 
un does not converge to uo in L 2 (T d ). That is, there exist e > and a subsequence 
{uN k ) such that 

\\uN k — uo||L2( T d) > e, 
for all fc. By Rellich-Kondrachov Theorem (Proposition 12.9ft . we have that there 
exists v G L 2 (T d ) and a further subsequence (also denoted by UN k ) such that 

u Nk k ^v, ini 2 (T d ). 

This implies that 

UN k — > W| weakly in L 2 (T d ), 
but this is a contradiction, since 

ujv fc — > uo, weakly in L 2 (T d ), 
and || 17 — uollL 2 (T d ) ^ e - Therefore, un — ► w in L 2 (T d ). The proof thus follows from 



expression (|34j) . □ 

This Proposition shows that even though the H- convergence only requires weak 
convergence in its definition, it yields a convergence in the strong sense (convergence 
in the L 2 -norm). 

6.2. Random environment. In this subsection we introduce the statistically ho- 
mogeneous rapidly oscillating coefficients that will be used to define the random 
VF-generalizcd difference elliptic operators, where the VY-generalized difference el- 
liptic operator was given in Section [5j 

Let (Q, fx) be a standard probability space and {T x : fi — ► D,;x £ Z d } be a 
group of ^-measurable and ergodic transformations which preserve the measure /x: 

• T x : Q, — > SI is ^"-measurable for all x £ Z d , 

• fJ,(T x A) = /i(A), for any A £ T and x £ Z d , 

• To = I , T x o T y = T x+y , 
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• For any / £ L 1 (f2) such that f(T x ui) = f(ui) /Lt-a.s for each x £ 7L d , is equal 
to a constant /x-a.s. 

The last condition implies that the group T x is crgodic. 

Let us now introduce the vector- valued ^-measurable functions {dj(u>);j = 
1, . . . , d} such that there exists 9 > with 

IT 1 < aj{w) < 9, 

for all uj £ O and j = 1, . . . ,d. Then, define the diagonal matrices A N whose 
elements are given by 

(35) a jj( x ) '■= a f = a 3 (T Nx uj) , x £ T N , j = 1, . . . , d. 

6.3. Homogenization of random operators. Let A > 0, Jn be a functional on 
the space of functions Hn ■ — > R, / £ H^~(T d ) (see also, subsection 12. 5|) . um 
be the unique weak solution of 

\u N - V N A N V^u N = f N , 

and uo be the unique weak solution of 

(36) Au - VAV w u = f. 

For more details on existence and uniqueness of such solutions see Sections [3] and 

m 

We say that the diagonal matrix A is a homogenization of the sequence of random 
matrices A N if the following conditions hold: 

• For each sequence Jn — > / in (T d ), converges weakly in H\,w to 
uo, when N — > oo; 

• a^dyy.u N — > 0^9^^, weakly in -k^^. (T d ) when N — > oo. 

Note that homogenization is a particular case of iJ-convergence. 
We will now state and prove the main result of this Section. 

Theorem 6.2. Let A N be a sequence of ergodic random matrices, such as the 
one that defines our random environment. Then, almost surely, A N (oj) admits a 
homogenization, where the homogenized matrix A does not depend on the realization 
u. 

Proof. Fix / £ i/~ 1 (T d ), and consider the problem 

Xu N - V N A N \7^u N = f. 

Using Lemma 15.31 and Remark 15.51 there exists a unique weak solution un of the 
problem above, such that its H^ w norm is uniformly bounded in N. That is, there 
exists a constant C > such that 

\\ u n\\h 1iW (t%) < C|l/llL 2 (T d )- 

Thus, the I/ 2 (T^)-norm of af d^.u^ is uniformly bounded. 

From VF-interpolation (see subsection 15. ip and the fact that Hi t w0^ d ) i s a 
Hilbcrt space (Lemma 12. 4|) , there exists a convergent subsequence of un (which 
we will also denote by wjv) such that 

Un — * uq, weakly in H\ t w (T d ), 
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and 

(37) afd^.u N ^VQ weakly in L 2 {T d ), 

as N — > oo; v being some function in L 2 x y^ w . (T d ). 

First, observe that the weak convergence in Hi ty(T d ) implies that 

(38) dw i UN Nj ^d Wi u weakly in L x ^ w .(T d ). 

From Birkhoff's ergodic theorem, we obtain the almost sure convergence, as N 
tends to infinity, of the random coefficients: 

(39) of—* Oi, 

where a.; = J3[af°], for any N Q e N. 

From convergences in (p?T|) . (f3"5)) and (j3"5)) . we obtain that 

v = aidwiUo, 

where, from the weak convergences, uo clearly solves problem (1351) . 

To conclude the proof it remains to be shown that we can pass from the sub- 
sequence to the sequence. This follows from uniquenesses of weak solutions of the 
problem ((36]). □ 



Remark 6.3. At first sight, one may think that we are dealing with a very special 
class of matrices A (diagonal matrices). Nevertheless, the random environment for 
random walks proposed in |14[ Section 2.3], which is also exactly the same random 
environment employed in (7j ; results in diagonal matrices. This is essentially due 
to the fact that in symmetric nearest-neighbor interacting particle systems (for ex- 
ample, the zero-range dynamics considered in a particle at a site x £ T N may 
jump to the sites x ± ej, j = 1, . . . , d. In such a case, the jump rate from x to x + ej 
determines the jth element of the diagonal matrix. 

Remark 6.4. Note that if u € Hw is a strong solution (or weak, in view of Remark 

El of 

Am - VAV w u = f 
and un is strong solution of the discrete problem 

Xu N - X7 N A N V^u N = f 

then, the homogenization theorem also holds, that is, un also converges weakly in 
Hi t w to u. 

7. Hydrodynamic limit of gradient processes with conductances in 

random environment 

Lastly, as an application of all the theory developed in the previous sections, 
we prove a hydrodynamic limit for gradient processes with conductances in random 
environments. Hydrodynamic limits for gradient processes with conductances have 
been obtained in for the one-dimensional setup and in |16j for the d-dimcnsional 
setup. However, the proof given here is much simpler and more natural, in view of 
the theory developed here, than the proofs given in [51ll6j. Furthermore, the proof 
of this hydrodynamic limit also provides an existence theorem for the VF-generalizcd 
parabolic equations in (]22p. 

The hydrodynamic limit allows one to deduce the macroscopic behavior of the 
system from the microscopic interaction among particles. Moreover, this approach 
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justifies rigorously a method often used by physicists to establish the partial differ- 
ential equations that describe the evolution of the thermodynamic characteristics 
of a fluid. 

This Section is structured as follows: in subsection 7.1 we present the model, 
derive some properties and fix the notations; subsection 7.2 deals with the hydro- 
dynamic equation; finally, subsections 7.3 and 7.4 are devoted to the proof of the 
hydrodynamic limit. 

7.1. The exclusion processes with conductances in random environments. 

Fix a typical realization u> £ f2 of the random environment defined in Section [6] 
For each x £ and j = 1, . . . , d, define the symmetric rate £ x x + e . — £, x +e x by 
(40) 

= a[(x) 

We, N \w((x + ej )/N)-W(x/N)] NlWjifa + iyty-Wifa/N)]' 

where (x) is given by (|33|) . and e\, . . . , e& is the canonical basis of M. d . Also, let 
b > -1/2 and 

Cx.x+ejiv) = 1 + b{q(x - e } ) + r)(x + 2 e,)} , 

where all sums are modulo N . 

Distribute particles on in such a way that each site of Tfj is occupied at 
most by one particle. Denote by 77 the configurations of the state space {0, 1} T « so 
that rj(x) =0 if site x is vacant, and ij(x) = 1 if site x is occupied. 

The exclusion process with conductances in a random environment is a 
continuous-time Markov process {r/ t : t > 0} with state space {0, 1} T « = {j] : 
— > {0, 1}}, whose generator acts on functions / : {0, 1} T « — > R as 

d 

(41) (L N f)( V ) = E E ^ x +e j c x , x+ej {r ] ) {f(a x ^ V ) - f( V )} , 

where a x ' x+£j r) is the configuration obtained from 77 by exchanging the variables 
r](x) and r/(x + e 3 -): 

{r)(x + ej) iiy = x, 
V (x) ify = x + e j , 

r](y) otherwise. 

We consider the Markov process {rj t : t > 0} on the configurations {0,1} T ™ 
associated to the generator in the diffusive scale, i.e., Lm is speeded up by TV 2 . 

We now describe the stochastic evolution of the process. After a time given 
by an exponential distribution, a random choice of a point x G is made. At 
rate Cx,K+ej the occupation variables 77(0;), r/(x + ej) are exchanged. Note that only 
nearest neighbor jumps are allowed. The conductances are given by the function W, 
whereas the random environment is given by the matrix A N := (a^(x)) dxd . The 
discontinuity points of W may, for instance, model a membrane which obstructs 
the passage of particles in a fluid. For more details see [16]. 

The effect of the factor c XiX + ej (rj) is the following: if the parameter b is positive, 
the presence of particles in the neighboring sites of the bond {x,x + ej} speeds up 
the exchange rate by a factor of order one, and if the parameter b is negative, the 
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presence of particles in the neighboring sites slows down the exchange rate also by 
a factor of order one. More details are given in Remark 17. II below. 

The dynamics informally presented describes a Markov evolution. A computa- 
tion shows that the Bernoulli product measures {v^ : < a < 1} are invariant, in 
fact reversible, for the dynamics. The measure is obtained by placing a particle 
at each site, independently from the other sites, with probability a. Thus, is a 
product measure over {0, 1} T « with marginals given by 

v£{r) : T}(x) = 1} = a 

for x in T d N . For more details see [H chapter 2]. 

Consider the random walk {X t }t>o of a particle in Tfj induced by the generator 
I/jv given as follows. Let £ x ,x+e.j given by ((40|) . If the particle is on a site x £ T%, it 
will jump to x + ej with rate N 2 ^ XjX+ej . Furthermore, only nearest neighbor jumps 
are allowed. The generator Ljv of the random walk {X t }t>o acts on functions 
/ : T% ->1 as 

=XX/ (£), 

3=1 

where, 

lw(^) =^{ ? „ + ,[/(£±a) -/(^)] +^,,[,(£^2) _,(*)]} 

It is not difficult to see that the following equality holds: 

d 

(43) h N f(x/N)=Y, d *M d wJ)(x) ■= V N A N ^f{x). 

3=1 

Note that several properties of the above operator have been obtained in Section 
[3] The counting measure on N~ 1 Tfj is reversible for this process. This random 
walk plays an important role in the proof of the hydrodynamic limit of the process 
rjt, as we will see in subsection 7.3. 

Let D(R + ,{0, 1} T «) be the path space of cadlag trajectories with values in 
{0, 1} T «. For a measure hn on {0, 1} Tjv , denote by P MJV the probability mea- 
sure on Z)(R_|_, {0, 1} T «) induced by the initial state and the Markov process 
{Vt '■ t > 0}. Expectation with respect to P MN is denoted by E MJV . 

Remark 7.1. The specific form of the rates c x ,x+ei * s n ot important, but two 
conditions must be fulfilled. The rates must be strictly positive, they may not depend 
on the occupation variables T)(x), T)(x + ej), but they have to be chosen in such a 
way that the resulting process is gradient, (cf. Chapter 7ifi| for the definition 
of gradient processes) . 

We may define rates c XjX + ei to obtain any polynomial $ of the form <£>(a) = 
a + Y^2<j<m a i a3 > m — 1> with 1 + E2<j<mi a i > 0- Let, for instance, m = 3. 
Then the rates 

Cx,x+ei{v) = Cx,x+e z (ri) + 

b{r](x - 2ei)n(x - e») + rj(x - ei)r\{x + 2e,) + r\{x + 2ei)rj(x + 3e,)} , 
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satisfy the above three conditions, where c x ,x+ei is the rate defined at the beginning 
of Section 2 and a, b are such that 1 + 2a + 36 > 0. An elementary computation 
shows that $(a) = 1 + aa 2 + ba 3 . 

7.2. The hydrodynamic equation. The hydrodynamic equation is, roughly, a 
PDE that describes the time evolution of the thermodynamical quantities of the 
model in a fluid. 

Let A = (ajj)dxd be a diagonal matrix with ajj > 0,j = l,...,d, and consider 
the operator 

d 

V.lVu : Y <h\ 

j=i 

defined on H>w- 

A sequence of probability measures {p N : N > 1} on {0, 1} T « is said to be 
associated to a profile po : T d — > [0, 1] if 



(44) \im^ N | E H(x/N)n(x) - J H(u)p (u)di 



>s y = o 



for every 5 > and every function H £ Bw 

Let 7 : T d — > [l,r] be a bounded density profile and consider the parabolic 
differential equation 



(45) 



dtp = VAVwHp) 
P(0,-) = 7(0 



where the function $ : [l,r] — > K is given as in the beginning of Section |U and 
t e [0,T], for T > fixed. 

Recall, from Section 01 that a bounded function p : [0,T] x T d — > [Z,r] is said 
to be a weak solution of the parabolic differential equation (|45[) if the following 
conditions hold. $(/£>(•,•)) and •) belong to L 2 ([0, T], ff w (T d )), and we have 
the integral identity 

p(t,u)H(u)du- [ p(0,u)H(u)du = f f $(p(s,u))VAV w H(u)duds , 

for every function H £ Hw and all t £ [0, T]. 

Existence of such weak solutions follow from the tightness of the process proved in 
subsection 7.3, and from the energy estimate obtained in Lemma l7.5l Uniquenesses 
of weak solutions was proved in subsection 14. II 



Theorem 7.2. Fix a continuous initial profile po : T d — > [0, 1] and consider a 
sequence of probability measures pn on {0,1} Tjv associated to po, in the sense of 
dill) . Then, for any t > 0, 



lim P MJV 

N— >oo 



Yl H{x/N) m {x)- f H(u)p{t,u)du 



>S \ = 



for every S > and every function H £ U>w Here, p is the unique weak solution 
of the non-linear equation (|45[) with I = 0, r = 1, 7 = po and $(a) = a + aa 2 . 
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Let M. be the space of positive measures on T with total mass bounded by one 
endowed with the weak topology. Recall that tt^ £ M. stands for the empirical 
measure at time t. This is the measure on T d obtained by rescaling space by N 
and by assigning mass 1/N d to each particle: 

(46) »f = pE th(x)6 x/n , 

where 5 U is the Dirac measure concentrated on u. 

For a function H : T d — > R, (ir^ , H) stands for the integral of H with respect to 

^,H) = ± H(x/N)r, t (x) . 

xef d N 

This notation is not to be mistaken with the inner product in L 2 (T d ) introduced 
earlier. Also, when ir t has a density p, w(t,du) = p(t,u)du. 

Fix T > and let D([0,T],Ad) be the space of .M-valued cadlag trajectories 
7r : [0, T] — ► M. endowed with the uniform topology. For each probability measure 
Hn on {0,1} T «, denote by Qjf the measure on the path space D([0,T], Ai) 
induced by the measure p n and the process tt^ 1 introduced in (J46[) . 

Fix a continuous profile po : T d — > [0, 1] and consider a sequence {pn ■ N > 1} 
of measures on {0, 1} T « associated to po in the sense (|44|). Further, we denote 
by Qw the probability measure on D([0,T], Ai) concentrated on the deterministic 
path n(t, du) = p(t, u)du, where p is the unique weak solution of (|45[) with 7 = po, 
Ik = 0, ffe = 1, k = 1, . . . , d and Q(a) = a + bo? . 

In subsection 17.31 we show that the sequence {Q^ N '■ N > 1} is tight, and in 
subsection 17.41 we characterize the limit points of this sequence. 

7.3. Tightness. The goal of this subsection is to prove tightness of sequence 
{^f^N N '■ N > 1}. We will do it by showing that the set of equicontinuous paths of 
the empirical measures (j46]) has probability close to one. 

Fix A > and consider, initially, the auxiliary .M-valued Markov process ' : 
t > 0} defined by 

n}< N (H) = (k? = ^ £ H»( X /N) Vt (x), 

for H in IS>w, where is the unique weak solution in H 1) w(^'n) (see Section [5]) 
of 

XHx - V N A N \7%-H^ = \H- VAV W H, 

with the right-hand side being understood as the restriction of the function to the 
lattice Tff (see subsection 15. 1[) . 

We first prove tightness of the process {n^ ,Af : < t < T},then we show that 
{n t A,JV : < t < T} and {n? : < t < T} are not far apart. 

It is well known [8] that to prove tightness of {n A,iv : < t < T} it is enough 
to show tightness of the real-valued processes {II t ' v (H) : < t < T} for a set of 
smooth functions H : T d — > R dense in C(T d ) for the uniform topology. 
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Fix a smooth function H : T d -> R. Keep in mind that U$' N (H) = (ir?,Hg), 
and denote by M^' X the martingale defined by 

(47) Mf< A = n^ N (H) - n x Q N (H) f : dsN*L N (n?,H») . 

Jo 

Clearly, tightness of IL^' (H) follows from tightness of the martingale M t ' and 
tightness of the additive functional J* ds N 2 Ln(k{? , 

A long computation, albeit simple, shows that the quadratic variation (M N ' X ) t 
of the martingale M^' X is given by: 



S-T E E fe<(^/iV)] 2 [^((x + e,-)/iV) - W^/JV)] x 

x / Cx >x+ej (rj s ) {ri s (x + ej) - rj s (x)} 2 ds 



3=1 x6T d 



In particular, by Lemma 15.31 

<^ A > t < ^Eii<ii 2 ^ < ^ 

3=1 

for some finite constant C{H), which depends only on H. Thus, by Doob inequality, 
for every A > 0, 5 > 0, 



(48) lim P„„ 

N— >oc 



sup \M t \ > 

0<t<T 



= 



In particular, the sequence of martingales {M^' X : N > 1} is tight for the uniform 
topology. 

It remains to be examined the additive functional of the decomposition (|47p. 
The generator of the exclusion process Ljy can be decomposed in terms of the 
generator of the random walk L/y. By a long but simple computation, we obtain 
that N 2 L n (tt n ,H^ ) is equal to 

Ei^ T,( h N^)(x/N) V (x) 

3 = 1 xET d N 

where {t x : x £ 1 d } is the group of translations, so that {r x rj)(y) = r){x + y) for 
x, y in Z d , and the sum is understood modulo N. Also, hij, hi,j are the cylinder 
functions 

Since is the weak solution of the discrete equation, we have by Remark 15.41 
that it is also a strong solution. Then, we may replace h^H^ by = XH^ — H 
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:s:s 



in the previous formula. In particular, for all < s < t < T, 



drN 2 L N (Tr?,H?) 



< 



(l + 3\b\)(t-s) 
N d 



It follows from the estimate given in Lemma [5.31 and from Schwartz inequality, that 
the right hand side of the previous expression is bounded above by C(H, b)(t — s) 
uniformly in N, where C{H, b) is a finite constant depending only on b and H. This 
proves that the additive part of the decomposition (|47|) is tight for the uniform 
topology and, therefore, that the sequence of processes {IT t ' : N > 1} is tight. 

N > 1} is tight for the uniform 



Lemma 7.3. The sequence of measures 
topology. 

Proof. Fix A > 0. It is enough to show that for every function H £ ]S>w and every 
e > 0, we have 



lim ] 

N — >oc 



sup \IL$> N (H) 

0<t<T 



n?,H)\>e 



0. 



whence tightness of ir^ follows from tightness of Ilf ' N . 
the last expression is bounded above by 



By Chcbyshcv's inequality, 



sup \TJi>»(H)-(-x?,H)\> 

0<t<T 



<2\\H?-H\\%, 



since there exists at most one particle per site. By Theorem 
O \\Hg - H\\% -> as N -> oo, and the proof follows. 



and Proposition 
□ 



7.4. Uniqueness of limit points. We prove in this subsection that all limit points 
Q* of the sequence Q^ N are concentrated on absolutely continuous trajectories 
7r(t, du) = p(t,u)du, whose density it) is a weak solution of the hydrodynamic 
equation (|4"5|) with I = 0, r = 1 and $(a) = a + ao? . 

We now state a result necessary to prove the uniqueness of limit points. Let, for 
a local function g : {0, 1} Z — ► R, g : [0, 1] — ► K be the expected value of g under 
the stationary states: 

g(a) = E, ya [ g (v)} ■ 

For I > 1 and <i-dimcnsional integer x = (xi, . . . , x<j), denote by rf{x) the empirical 
density of particles in the box B+(x) = {(yi, ■ ■ . , yd) £ Z d ; < y% — x% < £}: 

Proposition 7.4 (Replacement lemma). Fix a cylinder function g and a sequence 
of functions {F N : N > I}, F N : N' 1 ^ 



N 



such that 



lim sup -^j E F N{x/Nf < oo 



AT— too 



xETt 



Then, for any t > and any sequence of probability measures {/ijv : N > 1} on 

{o,i}K, 



lim sup lim sup E p 



N- 



w J2 F N (x/N){r x g( Vs )-~g( V f(x)) d s } 



= . 
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be a limit point of the sequence Q^ N and assume, without loss of 



The proof can be found in [16j subsection 5.3] 
Let 

generality, that Q^ N converges to Q* . 

Since there is at most one particle per site, it is clear that Q* is concentrated on 
trajectories ir t (du) which are absolutely continuous with respect to the Lebesgue 
measure, ir t (du) = p(t, u)du, and whose density p is non-negative and bounded by 
1. 



Fix a function H E H>w and A > 0. Recall the definition of the martingale M i 



N,\ 



introduced in the previous section. From 



we have, for every 5 > 0, 



lim 

N— >oc 



sup M t 

,0<t<T 



N,X\ 



> S 



= 



and from (|47|) . for fixed < t < T and 5 > 0, we have 

ft 



lim 



dsN*L N (-K«,H») 



> 5 



= 0. 



Note that the expression N 2 Ln(it^ , H^) has been computed in the previous 
subsection in terms of generator Ljy. On the other hand, L/v-ffj^ = XH^ — XH + 
\7AVwH. Since there is at most one particle per site, we may apply Theorem 
16.21 to replace {ir^,H^) and {ttq,H^) by (n t ,H) and (n ,H), respectively, and 
replace hjyH^ by VAVyyH plus a term that vanishes as N — > oo. 

Since E„ a [hij] = a 2 , i = 1, 2 and j = 1, . . . , d, we have by Proposition 17. 41 that, 
for every t > 0, A > 0, S > 0, i = 1, 2, 



lim lim sup 

e^0 N^oo 



Since r] e s (x) = e 
considerations, that 



1 ds W E VAV w H{x/N) 



X 



> 6 



0. 



^^dliLi [ X j/N, Xj/N + eej]), we obtain, from the previous 



lim lim sup < 



N — >oc 



(7T0 



> 5 



= 



Using the fact that Q]^' N converges in the uniform topology to Q*, we have that 



lim Q * 



(ir u G x H) - (ir ,G x H) 



f ds U(s- d n s (f[[- r 

Jo x i=1 



;])) . ^ 



> 5 



= o 



Recall that Q* is concentrated on absolutely continuous paths TTt(du) = p{t 1 u)du 
with positive density bounded by 1. Therefore, e~ d TT s (Ilj=i ['j ' Jr£e j\) converges in 
L x (T d ) to p(s, .) as e j 0. Thus, 



(7T t) ff) - (TT ,H) 



I ds m Ps ) , VAV W H) 
Jo 



> 5 



0. 
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.55 







sup < 


I 







Letting S I 0, we see that, Q* a.s., 

p(t,u)H(u)du- [ p(0,u)H(u)du = [ [ <$>(p(s,u))VAV w H(u)duds. 

T d Jf d Jo Jl d 

This identity can be extended to a countable set of times t. Taking this set to be 
dense we obtain, by continuity of the trajectories irt , that it holds for all < t < T. 

We now have a lemma regarding the energy of such limit points whose proof can 
be found in [T5], Section 6]. 

Lemma 7.5. There exists a finite constant K\, depending only on b, such that 

ds / dx (d X:) dw j H)(s 1 x)^(p{s,x)) 

-Kxf ds [ [dw^is^x^dix 1 ® Wj)\ < K . 
Jo Jr d J . 

From Lemma 17.51 we may conclude that all limit points have, almost surely, 
finite energy, and therefore, by LemmaO •)) S L 2 ([0, T], H hW (T d )). Anal- 

ogously, it is possible to show that /?(•,•) has finite energy and hence it belongs to 
L 2 ([0,T],iT w (T d )). 

Proposition 7.6. As N | oo, the sequence of probability measures converges 
in the uniform topology to Qw 

Proof. In the previous subsection, we showed that the sequence of probability mea- 
sures Q^ N is tight for the uniform topology. Moreover, we just proved that all limit 
points of this sequence are concentrated on weak solutions of the parabolic equa- 
tion (|45[) . The proposition now follows from the uniqueness proved in subsection 

o □ 



Proof of Theorem \ 7. S\ Since converges in the uniform topology to <Qw, a 

measure which is concentrated on a deterministic path, for each < t < T 
and each continuous function H : T d — > R, (tt^,H) converges in probability to 
J Td dup(t,u)H(u), where p is the unique weak solution of (|45[) with Ik = 0, r^ = 1, 
7 = po and $>(a) = a + aa 2 . □ 
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